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ABSTRACT 

Context. Near-Earth asteroid 162173 (1999 JU3) is the primary target of the Japanese Aerospace Exploration Agency (JAXA) 
Hayabusa-2 sample return mission, and is also on the list of potential targets for the European Space Agency (ESA) Marco Polo 
sample return mission. Earth-based studies of this object are fundamental to these missions. 
Aims. Our aim is to provide new constraints on the surface properties of this asteroid. 

Methods. We present a mid-infrared spectrum (5-38 yum) obtained with NASA's Spitzer Space Telescope in May 2008 and results 
from the application of thermal models. 

Results. These observations place new constraints on the surface properties of this asteroid. To fit our spectrum we used the near-Earth 
asteroid thermal model (NEATM) and the more complex thermophysical model (TPM). However, the position of the spin-pole, which 
is uncertain, is a crucial input parameter for constraining the thermal inertia with the TPM; hence, we consider two pole orientations. 
First is the extreme case of an equatorial retrograde geometry from which we derive a rigorous lower limit to the thermal inertia of 
150 Jm~ 2 s~°- 5 K~ [ . Second, when we adopt the pole orientation of Abe et al. (2008a) our best-fit thermal model yields a value for 
the thermal inertia of 700 ± 200 Jm~ 2 s~~ 05 K~ l and even higher values are allowed by the uncertainty in the spectral shape due to the 
absolute flux calibration. Our best estimates of the diameter (0.90 ± 0.14 km) and geometric albedo (0.07 ± 0.01) of asteroid 162173 
are consistent with values based on previous mid-infrared observations. 

Conclusions. We establish a rigorous lower limit to the thermal inertia, which is unlikely but possible, and would be consistent with 
a fine regolith similar to wthat is found for asteroid 433 Eros. However, the thermal inertia is expected to be higher, possibly similar 
to or greater than that on asteroid 25143 Itokawa. An Accurately determining the spin-pole of asteroid 162173 will narrow the range 
of possible values for its thermal inertia. 
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1. Introduction 

Near-Earth asteroid 162173 (1999 JU3) is the primary target of 
the Japanese Aerospace Exploration Agency (JAXA) Hayabusa- 
2 sample return mission. This object is also on the list of poten- 
tial targets for the European Space Agency (ESA) Marco Polo 
mission, which was selected for an assessment phase study in 
late 2007. The goal of both missions is to return to Earth sam- 
ples of a primitive near-Earth asteroid for laboratory analysis. 

Since these space missions seek to understand the origin and 
nature of volatile and organic material in the early Solar System, 
samples of a primitive asteroid (C, P, or D asteroids in the Tholen 
classification system; Tholen and Barucci 1989) are highly de- 
sirable. Asteroid 162173 is a C-class object, more specifically a 
Cg-type in the SMASS classification system (Binzel et al. 2001 ). 



In the asteroid population, C-class objects are abundant in the 
middle to outer Main Belt, and are primitive, volatile-rich rem- 
nants of the early solar system. Spectroscopy of C-class objects 
and comparisons with potential meteorite analogs suggest sur- 
face constituents of anhydrous silicates, hydrated minerals, or- 
ganic polymers, magnetite, and sulfides (e.g., Rivkin et al. 2002; 
Gaffey et. al.|2002). Dynamically, asteroid 162173 is most likely 
to have come from the main asteroid belt through the 3:1 jo- 
vian mean motion resonance (e.g., Bottke et al. 2002]!. At visi- 
ble and near infrared wavelengths (0.5 to 2.0 //m), this asteroid 
shows a spectrum with a neutral to bluish spectral slope (Binzel 
et al. 120011 Abe et al., |2008bb . A broad and shallow absorp- 
tion, centered near 0.7 fim and associated with hydrated minerals 
(iron-bearing phyllosilicates), has been reported at one rotational 
phase by Vilas (2008). This asteroid appears to be nearly spher- 
ical with a rotational amplitude of about 0. 1 magnitudes (Abe et 
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al. 2008b) and a rotation period of 0.3178 days (Hasegawa et al. 
2008). Ground-based and space-based narrowband mid-IR pho- 
tometry were used by Hasegawa et al. (2008 ) to derive estimates 
of the geometric albedo (0.063 +0.020 -0.015), diameter (0.91 + 
0.12 km), and thermal inertia (> 500 Jm^s^K -1 ). The small 
rotational amplitude makes determining a spin-pole orientation 
for this asteroid a challenging task; nevertheless, a spin-pole ori- 
entation has been estimated at ecliptic longitude and latitude of 
331 and 20°, respectively (Abe et al. 2008a). In this work, we 
present new mid-infrared spectroscopy (5-38 /im) obtained with 
NASA's Spitzer Space Telescope, which has allowed us to place 
new constraints on the surface properties of this asteroid. 

2. Observations and data reduction 

The spectrum was obtained with the Infrared Spectrograph 
(IRS) instrument (Houck et al. 120041 ) on NASA's Spitzer Space 
Telescope (Werner et al. 120041 . on UT May 2.084, 2008. The 
asteroid was observed at a heliocentric distance of 1.202 AU, 
at a Spitzer distance of 0.416 AU, a phase angle of 52.6°, and 
a galactic latitude of 3.2°. The IRS measures spectra over the 
range 5.2-38 ^im. The observations presented here used the low 
spectral resolution mode (R = A/AA ~ 64-128), which covers 
the full spectral range in four long-slit segments, and the final 
spectrum corresponds to a total on-object exposure time of 738.2 
seconds. These four segments are short wavelength, 2nd order 
(SL2; 5.2-8.5 yum), short wavelength, 1st order (SL1; 7.4-14.2 
/mi), long wavelength, 2nd order (LL2; 14.0-21.5 fim), and long 
wavelength, 1st order (LL1; 19.5-38.0 pm). When each of the 
low-resolution orders is centered on the target, the other is offset 
(~ 50" for SL and ~108" for LL) from the target, providing a 
measurement of the sky background. We anticipated a crowded 
field, so we also obtained the sky background from a "shadow" 
observation, which duplicated the first observation but with the 
asteroid out of the field of view. We did not use the peak-up rou- 
tine because (a) the ephemeris uncertainty was only about 0.3 
arcsec, much less than the slit widths, (b) the Spitzer pointing 
accuracy is about 1 arcsec, also less than the slit widths, and (c) 
the object's Galactic latitude was low so we did not want IRS to 
inadvertently peak-up on a background source. Our spectra were 
cleaned, extracted, and co-added using the Spitzer IRS Custom 
Extraction tool (SPICE), and our own routines. More details on 
our data reduction can be found in Kelley et al. (2009). 

When the reduction was completed, we found a discrepancy 
of 10% between the fluxes of overlapping wavelengths in the SL 
and LL orders, which is equal to the expected uncertainty in the 
absolute flux calibration (Houck et al. 2004), and small offsets 
like this one are not unusuafj We scaled the LL and SL fluxes 
to match each other in the overlapping wavelengths, the SL flux 
by 1.05 and the LL by 0.95, and we include this 10% system- 
atic absolute calibration uncertainty in the modeled parameters 
discussed in the next section. This spectrum is plotted in Fig. 1 
along with the best-fit thermophysical model (TPM, Sect. 3). 

3. Thermal modeling 

The flux from near-Earth asteroids in the wavelength range cov- 
ered by the IRS is dominated by thermal emission. The mea- 
sured spectral energy distribution (SED) depends on the object's 
size, composition, and temperature distribution. This last term is 

1 We considered and ruled out other potential sources for the SL and 
LL flux difference, such as placement of the object within the slit and 
rotational variability of the asteroid. 




Wavelength (iim) 

Fig. 1. The Spitzer spectrum of near-Earth asteroid 162173 
( 1 999 JU3) plotted along with our best-fit thermophysical model 
(TPM). The model parameters are given in Table 1 and discussed 
in Sect. 3. The LL and SL fluxes have been scaled to match each 
other in the overlapping wavelengths (see Sect. 2). 



dependent on several factors, including distance from the Sun, 
albedo, thermal inertia, surface roughness, rotation rate, shape, 
and spin-pole orientation. We used two models to fit our 5 to 38 
pm spectrum, the relatively simple near-Earth asteroid thermal 
model (NEATM, Harris 1998) and the more complex thermo- 
physical model (TPM, e.g., Harris and Lagerros 2002 and refer- 
ences therein). 

The NEATM is a refinement of the standard thermal model 
(STM, Lebofsky et al. [19861 Lebofsky and Spencer [19891 , 
which was developed and calibrated for main-belt asteroids. The 
NEATM accounts for observations at larger phase angles and 
contemplates possible deviations from the STM in the thermal 
behavior of smaller asteroids. Unlike the STM, the NEATM re- 
quires observations at multiple wavelengths and uses this infor- 
mation to force the model temperature distribution to be con- 
sistent with the apparent color temperature of the asteroid. The 
NEATM solves simultaneously for the beaming parameter (if) 
and the diameter (D). The beaming parameter was originally in- 
troduced in the STM to allow the model temperature distribution 
to fit the observed enhancement of thermal emission at small so- 
lar phase angles due to surface roughness. In practice, rj can be 
thought of as a modeling parameter that allows a first-order cor- 
rection for any effect that influences the observed surface tem- 
perature distribution (such as beaming, thermal inertia, and ro- 
tation). The values of rj derived for asteroid 162173 using the 
NEATM fits to our spectrum are given in Table 1; for all our 
thermal models we assumed a bolometric emissivity (e) of 0.9. 
These rj values can be compared with those of other asteroids 
studied in a similar manner, and in fact fit well the trend of in- 
creasing T] (lower thermal fluxes and cooler color temperatures) 
with increasing solar phase angle, observed for about 40 near- 
Earth asteroids so far (e.g., Delbd 2004, Campins et al. 2009). 
We did not make any observations of this asteroid's reflected 
continuum, so to estimate the geometric albedo (pv), we adopted 
the values given by Hasegawa et al. ( [20081 of H v = 18.81+0.03 
and G = -0.1 15+0.009, where Hy is the visible absolute magni- 
tude and G is the slope parameter. 

In Table 1, scaling means that a correction was applied to 
the fluxes in the LL and SL orders (Sect. 2). Diameter, jj, F, 
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Model 


Scaling 


Diameter 


n 


r 


Pv 






(km) 




(Jm- 2 s- - 5 K- 1 ) 




NEATM 


no scaling 


0.97 + 0.15 


1.90 + 0.17 




0.06 + 0.01 


NEATM 


scaled orders 


0.91 +0.14 


1.63+0.15 




0.07 ± 0.01 


TPM1 


scaled orders 


0.87+0.13 




150 


0.07 ± 0.01 


TPM2 


no scaling 


0.97 + 0.15 




-1500 (Fig. 2) 


0.06 + 0.01 


TPM2 


scaled orders 


0.90 + 0.14 




700 ± 200 


0.07 + 0.01 



and p v are the derived diameter, beaming parameter, thermal 
inertia, and geometric albedo, respectively. TPM1 yields a rig- 
orous lower limit to the thermal inertia (Sect. 3) and TPM2 as- 
sumes the pole orientation of Abe et al. (|2008a) given in Sect. 
1 . Uncertainties in the diameter estimated with thermal model- 
ing usually exceed the formal errors and are typically 10-15% 
(e.g., Delbd et al. 2003). The listed error in the geometric albedo 
includes the uncertainty in Hy, with a 0.05 magnitude compo- 
nent due to the rotational amplitude of 0. 1 magnitudes, and the 
diameter uncertainty propagated into the albedo. 

The TPM is applied to asteroids with known shape and spin 
state, and is generally used to constrain the thermal inertia (F), 
which is particularly relevant to the proposed sample-return mis- 
sions since F can be diagnostic of the regolith. For example, typ- 
ical values of F (in S.F units Jm~ 2 s~ 5 K~ l ) are 50 for lunar re- 
golith, 400 for coarse sands, and 2500 for bare solid rocks (e.g., 
Delbo and Tanga [20091 and references therein). As mentioned 
in Sect. 1, there are estimates of the rotation period, approxi- 
mate spin-pole orientation and the shape of asteroid 162173, so 
in principle we can apply the TPM. However, because of the un- 
certainties in the pole position, the TPM does not uniquely con- 
strain the thermal inertia of this asteroid; we elaborate on this 
point in the next paragraph. Another input parameter needed in 
this case is the surface roughness of the asteroid. This roughness 
can also be constrained by the TPM if observations are made at 
several phase angles, but in our case it has to be assumed be- 
cause we obtained only one spectrum. Little is known about the 
typical surface roughness of NEAs, and there is no consensus 
about what roughness parameters should be used (e.g., Harris et 
al 2007 and references therein). We explored the effect of surface 
roughness by incorporating two values of surface roughness into 
our models (see below). The difference in the derived thermal 
inertia between the smooth surface case and the rough case was 
less than 30%. Because this uncertainty is considerably smaller 
than that due to the spin-pole orientation, it is not meaningful to 
try to refine the roughness parameters at this time. 

As mentioned, the spin-pole orientation is a crucial input pa- 
rameter for constraining the thermal inertia. However, the esti- 
mate of the spin-pole position for asteroid 162173 is uncertain 
because the amplitude of the rotational light curve is small and 
the range of geometries for the light curve observations has been 
limited. We, therefore, consider two cases for the orientation of 
the spin-pole. The first is the extreme case of an equatorial ret- 
rograde geometry, i.e., the spin-pole pointing at the South eclip- 
tic pole. As input parameters for this TPM, we used a spherical 
shape (a reasonable assumption given the small amplitude of the 
light curve), the equatorial retrograde geometry, the rotation pe- 
riod of 0.3178 days, and a perfectly smooth surface (this last 
parameter yields a lower thermal inertia than in the rough sur- 
face case). Then we allowed the unconstrained parameters in the 
model (D, py, and F) to vary, an approach that worked particu- 
larly well because our spectrum spans a broad wavelength range 
that includes the peak in thermal emission. We generated a grid 



of models with F ranging from to 2500 Jm~ 2 s~°- K~ x , and the 
best-fit model was determined by minimizing the reduced^ 2 be- 
tween the model fluxes and our data. This case (equatorial ret- 
rograde geometry and smooth surface) yields a rigorous lower 
limit to the thermal inertia of 150 Jm~ 2 s~ 5 K~ l . 

For the second case, we adopted the pole orientation of Abe 
et al. (2008a) given in Sect. 1, which implies a subsolar latitude 
of approximately 30 at the time of the Spitzer observations. As 
in the first case, we generated a grid of models, and the best-fit 
model was determined by minimizing the reduced^ 2 . In Fig. 1, 
we plot our best-fit TPM model along with the observed spec- 
trum, and the corresponding parameters are listed in Table 1, 
again assuming e = 0.9 and the Hy and G from Hasegawa et al. 
(2008). Fig. 2 shows the reduced^ 2 vs thermal inertia and the es- 
timated diameter and geometric albedo vs thermal inertia for the 
rough-surface case (the parameters for the rough-surface model 
are y — 53° [half-angle opening of spherical section craters] and 
f = 0.77 [fractional surface coverage of craters]). The smooth 
surface model is not plotted because of space limitations, but fits 
equally well and yields about 30% lower values for the thermal 
inertia. 

Uncertainties in thermal modeling usually exceed the formal 
errors resulting from the scatter of the flux measurements. The 
main sources of uncertainties in our derived physical parame- 
ters are the absolute photometric uncertainty in the thermal flux 
(diameter), the relative spectral uncertainty (temperature, beam- 
ing parameter and thermal conductivity), and the uncertainty in 
Hy (albedo), including the rotational light curve. Comparison 
of diameters derived from thermal models with those derived 
from other sources, such as radar, indicates that overall errors 
are normally less than 15% in diameter (Delbo et al. [2003). 
The contribution of relative spectral uncertainty to errors in the 
thermal continuum and derived parameters is dominated by the 
absolute-flux-calibration adjustment of the LL and SL segments 
of the spectrum (Sect. 2). The estimated uncertainties are given 
in Table 1, and in the case of the (second case) TPM, they are 
also illustrated in Fig. 2. 

Spectral features superposed on the mid-infrared thermal 
continuum have been detected in other asteroids and can be com- 
positionally diagnostic. However, except for a possible weak fea- 
ture near 15 pm (Fig. 1), no obvious structure is apparent in our 
spectrum, and a detailed analysis is deferred to a later publica- 
tion. 



4. Results and conclusion 

Our results are summarized in Table 1, and the most relevant 
to the proposed sample return missions are the constraints our 
observations impose on the thermal inertia, F. When we adopt 
the extreme case of an equatorial retrograde geometry and a 
smooth surface for this asteroid, we derive a rigorous lower limit 
to the thermal inertia of 150 JmT 2 s~ a5 K~ l . This value would 
be consistent with a fine regolith, similar to the one found in 
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Fig. 2. The the reduced^ 2 vs thermal inertia for a rough surface 
TPM. The parameters for the rough surface model are given in 
Sect. 3. The bottom plot shows size and geometric albedo vs 
thermal inertia. The best-fit model was determined by minimiz- 
ing the reduced^ 2 between the model and the data, and is shown 
in Figure 1 along with the observed spectrum. 



asteroid 433 Eros. Adopting the pole orientation of Abe et al. 
(2008a) yields a best-fit model with thermal inertia at 700 +/- 



200 Jm 



-2 „-0.5 v-\ 



K . Even higher values are allowed by the un- 



certainty in the spectral shape due to absolute flux calibration. 
This best-fit thermal inertia and its size place asteroid 162173 
about two-sigma above the trend toward increasing thermal in- 
ertia with decreasing asteroid diameter as discussed by Delbo 
et al. ( 120071 1 . The range in possible thermal inertias for asteroid 
251 43 Itok awa adopted by Delbo et al. d2007l l (based on Muller 
et al.|2005]and MuellerEEHl is 350 to 800 Jm^s^K- 1 , which 
is within the range of possible values for asteroid 162173. If the 
spin-pole of Abe et al. d2Q08al) turns out to be correct, these last 
two results would suggest that the surface of asteroid 16217 has a 
texture rockier than or similar to that of asteroid 25143 Itokawa. 
However, given the uncertainties in the spin-pole position, we 
hesitate to even make this analogy. 

Our estimates of the diameter and albedo of asteroid 162173 
are consistent with those of Hasegawa et al. (2008 ). However, 
there are marginally significant discrepancies that may deserve 
further attention. There are temperature differences between 
their space-based (Akari) and ground-based (Subaru) data sets 
and between their data and ours. As a way to illustrate these dif- 
ferences we give our best-fit NEATM parameters for their Akari 
and Subaru data, which are D = 0.80+0. 12 km, p v = 0.08+0.03, 
77= 1 .0+0.4 and D = 1 . 1 3±0. 17 km, p v = 0.04+0.02, 7 = 2.1 +0.6, 
respectively. Only for r\ are the differences slightly larger than 



one sigma. It may be possible to constrain the spin-pole orien- 
tation better by modeling all three data sets (ours and the two 
sets in Hasegawa et al. 2008 ) simultaneously, allowing the pole 
solution to be a free parameter; however, this is work beyond the 
scope of this letter. 

In summary, thermal models of our Spitzer spectrum of as- 
teroid 162173 (1999 JU3) provide a rigorous lower limit to the 
thermal inertia of 150 Jm 1 s~ 5 K 1 , which is unlikely but pos- 
sible. This value would be consistent with a fine regolith, similar 
to the one found in asteroid 433 Eros. However, if we adopt the 
pole orientation of Abe et al. (2008a), the best-fit model yields a 
thermal inertia of 700 + 200 Jm~ 2 s~ 5 K~ l . Even higher values 
are allowed by the uncertainty in the spectral shape due to abso- 
lute flux calibration. Our best estimates of the diameter (0.90 ± 
0.14 km) and geometric albedo (0.07 + 0.02) of ast eroid 162173 
are consistent with those of Hasegawa et al. (2008 ). 
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